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Abstract. I present results on the interstellar medium towards the O stars observed in 
five optical spectroscopic surveys: GOSSS, OWN, lACOB, NoMaDS, and CAFE-BEANS. 
I have measured both the amount [£(4405 - 5495)] an d type [£ 5495 ] of extinction to wards 
several hundreds of Galactic O stars and verified that the lMaiz Apellaniz et alj ( l2014ah fam¬ 
ily of extincti on laws provides a si gnifica ntly better fit to op tical-rNIR Galactic extinction 
than either the lCardelli et alj ( Il989ll or the lFitznatrickl d 19991) families. £5495 values are con¬ 
centrated between 3.0 and 3.5 but for low values of £(4405 - 5495) there is a significant 
population with larger £5495 associated with Hit regions. I have also measured different 
DIBs and I have found that IV(5797)/W(5780) is anticorrelated with £5495, a sign that ex¬ 
treme f clouds are characterized not only by low ionization environments (as opposed to 
<T clouds) but also by having a larger fraction of small dust grains. The equivalent width 
of the “Gaia DIB” (8621 A) is strongly correlated with £(4405 - 5495), as expected, and 
its behavior appears to be more cr-like than ^-like. We have also started analyzing some 
individual sightlines in detail. 
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1. 0-Star spectroscopic surveys 


There are several ongoing ground-based 
optical spectroscopic surveys of early-type 
stars (mostly of O spectral type but also in¬ 
cluding B and WR stars) whose main goal is 
to study the properties of the stars. However, 
those surveys include large amounts of in¬ 
formation on the i ntervening ISM, which we 
are also analyzing (Penades Ordaz et al.ll201^ 


iMafz Apellaniz et al.ll2014b). 


The G alactic O-Star Spectroscopic Survey 
(GOSSS, iMaiz Apellaniz et ^ 1201 ih is ob¬ 
taining blue-violet, R ~ 2500, high S/N spec¬ 
troscopy of all O-star candidates that are ac¬ 


cessible to the survey telescopes. Results have 
already been pub l ished in two survey papers 
dSota et al.l 2()llL 2()14t) a nd elsewhere (e.g. 
IWalborn et al.ll2010l 1201 ll) . To date, we have 
observed and processed 2087 stars and several 
hundred more have been observed but are still 
unprocessed. The processed sample includes 
~900 O stars and ~900 B stars up to Ay ~ 12. 
The faintest stars are being observed with the 
11m GTC in the north and the 8 m Gemini in 
the south. 


GOSSS is reaching a large number of stars 
but it is limited by its intermediate resolution, 
coverage of only the blue-violet region (though 
some of the spectra reach longer wavelengths), 
























2 


Maiz Apellaniz: The ISM in O-star spectroscopic surveys 


[-•H Jo+2M 
h»H Jo+2M+St 
[-•H Jo+2M+Ty 
h«H Jo+2M+Ty+St 
O Oe 



0.0 0.5 1.0 1.5 2.0 2.5 

E(4405-bA9S) 


h*H Jo+2M 
hn Jo+2IVI+St 
[-•H Jo+2M+Ty 
hn Jo+2M+Ty+St 
O Oe 



0.0 0.5 1.0 1.5 2.0 2.5 

£(' 4405 - 5495 ) 


hn Jo+2M 
[-•H Jo+2M+St 
i-»i Jo+2M+Ty 
h«-i Jo+2M+Ty+St 
O Oe 



0.0 0.5 1.0 1.5 2.0 2.5 

£(' 4405 - 5495 ) 


Fig. 1. xii of the CHORIZOS fits as a function of £(4405 - 5495) for the GOSSS I+II sample using the 
F99 [top left], CCM [top right], and M14 [bottom] families of extinction laws. Note that the vertical scale 
is the same for the three plots. Different symbols are used for stars with dilferent sets of photometric data 
(Jo = Johnson, 2M = 2MASS, Ty = Tycho-2, St = Strdmgren). Oe stars (which are expected to have poorer 
fits due to IR excesses) and stars with £ 5495 > 4 are marked. 
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Fig. 2. i ?5495 as a function of £(4405 - 5495) for the 
CHORIZOS fits for the GOSSS I+II sample using 
the M14 family of extinction laws. See Fig. [T] for 
the symbol nomenclature. 


and single-or-few epochs available per star. 
There are several surveys that complement 
those limitations by obtaining high-resolution 
spectroscopy of a subsample of the GOSSS 
targets. OWN is observing the southern stars 
using three different telescopes and multi¬ 
ple epochs with the main purpose of study - 
ing spectroscopic binarity (iBarba et al.ll2010h . 
lACOB and lACOB-sweG are observing 
northern stars with the mai n purpose of char¬ 
acteri zing their properties dSimon-Dlaz et alJ 
I 2 OI ih . NoMaDS is extending lACOB towards 
fainter magnitude s using the 9 m Hobby- 
EbeiTy Telescope dMafz Anellaniz et al.ll2012l : 
iPellerin et al.l |2012|) . Finally, CAFE-BEANS 
is doing multi-epoch spectroscopy of north¬ 
ern stars to study their sp ectroscopic bina¬ 
rity dNegueruela et al.l[2015h . Combining those 
four high-resolution surveys, ~450 O stars and 
another ~450 early-type stars have been ob¬ 
served so far. 


2. Dust: amount and type 

We have recently derived a new family 
of optical-nNI R extinction laws based on 30 
Doradus data dMafz Anellaniz et al.ll2014^ . In 


that paper we used a small data set to check 
its validity for Galactic sightlines. Our first 
goal here is to extend that verification to the 
larger data set of 448 Galactic O stars of the 
first two GOSSS survey papers (GOSSS I-hII, 
ISota etalJlMl |20l4 . For that purpose we 
have used the Galactic O-Sto Catalog (GOS C, 
iMafz Anellaniz et ^l2004t ISota et alJl2008l to 

collect Johnson UBV and 2MASS JHK^ pho¬ 
tometry for all of the stars in the sample (af¬ 
ter merging some unresolved binary systems) 
and Tycho-2 BV and Stromgren uvby photome¬ 
try in the cases where they were available. The 
photometry was proce s sed usin g CHORIZOS 
dMafz ApellanizI [2004 l2013bh to obtain the 
amount [£(440 5 - 5495)1 and type [£5 4951 of 
extinction (see iMafz Anellanizl l2013al for an 
explanation of the choice of those quanti¬ 
ties) in_an_malogous way to that used by 
e.g.lArias et alJ d2006h or lMafz Anellaniz et alJ 
d2014al) . Note that this procedure can be car¬ 
ried out with high precision due to the use 
of the Teff derived from the GOSSS spectral 
types. 

The CHORIZOS experiments were exe¬ 
cuted three ti mes using the F99 dFitzpatrickl 
Il999h. CCM dCardelli etal.l ll989l) . and M14 
dMafz Apellaniz et alJl2014a) to test which of 
those families of extinction laws yields better 
fits to Galactic extinction. The values (ex¬ 
pected to be ~ 1 in the ideal case) are shown in 

Fig.m 


- The F99 results yield, on average, the poor¬ 
est fits. Although they are not bad for 
low values of £(4405 - 5495) and £ 5495 , 
they worsen when either of those quanti¬ 
ties increase (they are significantly bad for 
£5495 > 4). 

- The CCM results are better that those 
of F99 in terms of However, they 
yield poor results when Stromgren pho¬ 
tometry is used due to their use of a sev¬ 
enth degree polynomial for interpolating in 
wavelength (F99 and M 14 use splines, see 
iMafz ADellanizll2013ah . 

- The M14 results are the best of the three 
families. Indeed, the only case with;^^J^J > 
3 corresponds to an Oe star, which are 
expected to yield poor fits in some cases 
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Fig. 3. 1^5495 histograms for the (Galactic) GOSSS I+II [left] and (30 Doradus) VFTS [right] samples. 


due to the IR excesses produced by their 
disks. Therefore, the M14 family is the best 
choice for optical+NIR extinction not only 
in 30 Doradus but also in the Milky Way. 

After choosing the M14 family, I show its 
detailed CHORIZOS results in Figs. |2] and [3 
The latter also includes an R 5495 histogram 
for the VFT S O-star sa mple in 30 Doradus 
dEvans et al.l l201 ll : iMafz Anellaniz et~^ 
l2014ah . The main findings are: 

- As previously known for many decades, 
the majority of sightlines have R 5495 be¬ 
tween 3.0 and 3.5. 

- At very low extinctions (£(4405 - 5495) < 
0 . 2 ) the error bars on R 5495 are too large to 
yield significant results. 

- There are few stars with R 5495 <3.0 (sight¬ 
lines with a large proportion of small dust 
grains) but they are a significant fraction of 
the objects with large extinction. 

- In the range 0.2 < £(4405 - 5495) < 1.2 
there is a significant fraction of stars with 
large values of £ 5495 . These sightlines are 


depleted in small grains and are associated 
with Hii regions. 

- As previously noted, some Oe stars have a 
poor fit due to their IR excesses. 

- The £5495 histograms for the Galaxy and 
30 Doradus are markedly different, with 
the latter showing a larger fraction of 
high -£5495 sightlines. The differences can 
be explained by [a] the lower values of 
£(4405 - 5495) and [b] the larger fraction 
of Hit region sightlines in 30 Doradus. 

In the future we plan to extend the analysis 
to the rest of the GOSSS sample, which is on 
average more extinguished than the one in this 
work. 

3. DIBs 

Another line of research using these O-star 
surveys is the study of the Diffuse Interstellar 
Bands (DIBs). Our long-term plans include 
the analysis of several tens of DIBs but 
here I concentrate on three of them: the 
5780 and 5797 DIBs (central wavelengths of 
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Fig. 4. [left] IV(5780) as a function of £(4405 - 5495) and [right] W{5191) as a function of IT(5780) for 
the GOSSS sample. 


5780.48 A and 5797.06 A, respectively), which 
are among the most studied DIBs, and the 
8621 DIB (central wavelength of 8620.65 A), 
which is less studied but has received recent 
atention due to its inclusion in the wavelength 
region being observed by the Gaia RVS instru¬ 
ment. 

We have extracted the equivalent widths 
(W) for the three DIBs above using a subset of 
the high-resolution spectra in the surveys, cor¬ 
recting for saturation effects in the cases with 
deeper absorption profiles (this is not an issue 
for low values of extinction but some of our 
lines are heavily extinguished). These type of 
studies have been done before and the origi¬ 
nality of this ongoing work lies in: 

- The large sample in the survey, which in its 
final version will include over 1000 early- 
type star sightlines. 

- The large dynamic range in extinction, 
which will eventually extend to Ay ~ 12 
(note that the subsample presented here 
only extends to Ay ~ 6 ). 


- The completeness down to a given magni¬ 
tude (currently, Z? ~ 8 , to be extended in 
the future) for O stars and the full-sky cov¬ 
erage. 

- The use of modern values of 
£(4405 - 5495) and £5495 (see above) 
as extinction measurements. 

I plot in Fig. |4] VT(5780) as a function of 
£(4405 - 5495) and W{5797) as a function of 
1T(5780). In both cases the two quantities are 
highly correlated but the Pearson correlation 
coefficient (shown in both plots) is between 0.8 
and 0.9, a result similar to that of previous stud- 


ies (e.g. 

Friedm 

an et alJl 201 ll IVos et alJl 201 ll 

[Raimond et al.l 

2012). More interestinalv, in 


the right panel we can see that sightlines with 
low values of £5495 tend to lie above the av¬ 
erage relationship while those with high val¬ 
ues of £5495 tend to lie below the average rela¬ 
tionship. To better visualize the effect, I plot in 
Fig.| 5]£5495 as a function of VT(5797)/VT(5780). 
The Pearson correlation coefficient in that 
case is -0.384, significant but not very large. 
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Fig. 5. Ri 495 as a function of IV(5797)/IV(5780) for 
the GOSSS sample. 


However, the anticorrelation is clear when we 
examine the extremes of the distribution; all 
sightlines with 1F(5797)/VF(5780) > 0.50 have 
/?5495 < 3.5 and all with 7^5495 > 5.0 have 
1F(5797)/VF(5780) < 0.35. Putting it in an¬ 
other way, most sightlines concentrate in the 
lower left quadrant of Fig. |5] with some points 
in the lower right and upper left quadrants and 
none in the upper right. 

How do we interpre t this result? It has been 
known for some time dKrelowski et al.l [19971 : 
ICox et al.1 l2005h that DIBs come in families 
and that the 5780 and 5797 DIBs are represen¬ 
tative of two types of environments. The 5797 
DIB is (relatively) strong in ^ Oph, a sightline 
known to be exposed to a low UV flux, and is 
thus termed a DIB”. On the other hand, the 
5780 DIB is strong in cr Sco, a sightline ex¬ 
posed to a high UV flux, and is thus termed a 
“cr DIB”. See Fig. | 6 ] and notice the large dif¬ 
ference in the 5780 DIB for those two sight¬ 
lines, which have a similar £(4405 - 5495). 
( Oph and cr Sco are the two prototypical ex¬ 
amples of this phenomenon but in our sample 
we have found even more extreme cases such 
as HD 207 538 and 0 Mus B, also in Fig. |6l 
What we have found is that ^ clouds (those 
with low UV flux) tend to have low values 
of £5495 (i.e. a larger proportion of small dust 


grains) while cr clouds (those with high UV 
flux) tend to have large values of £5495 (i.e. a 
smaller proportion of small dust grains). Such a 
relationship between cloud typ e and £5495 was 
proposed bv ICami et al.l (Il997h . indicating that 
^ clouds likely correspond to cloud cores and cr 
clouds to cloud skins, but here we demonstrate 
its existence with a large sightline sample for 
the first time. 

The difference in grain population between 
denser, non-UV-exposed, ^ clouds and thinner, 
UV-exposed, cr clouds can be interpreted in 
terms of selective destruction of small grains 
in the latter. One interpretation would be that 
sputtering dominates over shattering in Hit 
regions and other parts of the ISM exposed 
to UV light, as the former preferentially de¬ 
stroys small grains and th e latter preferen¬ 
tially destroys large g rains dJones et al.lll994t 
[Andersen et akll^l lh. However, sputtering is 
only expected to become efficie nt above 10^ K 
for refractory grains (Fig. 25.4 in lDrainel201 ih . 
so a different mechanism is required. A possi¬ 
bility is that the extreme UV radiation present 
in Hit regions destroys the PAHs th at likely 
const itute the small grain population dTielensI 

I2nn8h . 

I also show in Fig. [T] VT(8621) as a func¬ 
tion of £(4405 - 5495). As expected, both 
quantities are strongly correlated with Pearson 
coefficients similar to those of other DIBs. 
Such measurements had been carried out be¬ 
fore for that DIB but reaching lower values of 
£(440 5 - 5495) dMiinari et al J2008 HKos et all 
l20I3h . In the future we will extend this study 
to even larger extinctions. 

4. Individual sightlines 

A third line of work is the study of indi¬ 
vidual sightlines in detail, of which we are an¬ 
alyzing several. Here we present a summary 
of our results on Berkeley 90, a little-studied 
cluster that includes two early-O-type systems, 
LS III -h 46 11 a nd LS III -h 46 12, separated b y 
a short distance dMafz Anellaniz et alJl20159) . 

Both stars show large extinctions but 
LS III -1-46 11 is significantly more extincted 
than LS III -1-46 12 (by 30%). £ 5495 , on the 
other hand is very similar for both stars and 
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Wavelength (A) Wavelength (A) 

Fig. 6. [left] The prototype ^ (f Oph) and cr (cr Sco) sightlines for the 5780 and 5797 DIBs and [right] two 
more extreme examples, HD 207 538 (() and 9 Mus B (cr). 



Fig. 7. 1T(8621) as a function of £(4405 - 5495) 
for the GOSSS sample. 

has standard values between 3.3 and 3.4. We 
have measured as many as seven atomic, seven 
molecular, and fifty DIBs for both stars and 


we have found large differences among them. 
We show in Fig. [8] the equivalent-width ra¬ 
tios for the lines we have measured. From the 
kinematic point of view, there are two com¬ 
ponents clearly visible in our high-resolution 
spectra. The weak one has a velocity interme¬ 
diate between that of the cluster and ours, ap¬ 
pears to be of relatively low density, and af¬ 
fects the two stars similarly. The strong one 
has a velocity similar to that of the clus¬ 
ter and affects the two stars differently: for 
LS III -1-46 12 it appears to be of lower den¬ 
sity than for LS HI -1-46 11. Combining all the 
information, we propose a model in which 
three different “clouds” (actually, one simple 
cloud and one cloud with core and skin) are lo¬ 
cated between Berkeley 90 and us (Fig. ©.The 
reader interested in more det ails is referred to 
iMafz Apellaniz et’aP (l2015ah . 

One interesting result of the Berkeley 90 
analysis is that it alows us to classify a large 
number of DIBs in a cr - f scale, something 
that was not possible before due to the lack of 
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Wavelength (A) 

Fig. 8. Equivalent-width ratios between 
LS III+46 12 and LS III+46 11. Red is used 
for atomic lines, green for molecular lines, and 
blue for DIBs. All points show the EW uncorrected 
for saturation (which should be a good approxi¬ 
mation in all cases except for the intense atomic 
lines) except for Ki 47664.911+7698.974. In that 
case, besides the uncorrected values we show the 
corrected one (as an unfilled symbol) calculated 
using a kinematic decomposition. The x value 
in each case is Tq, the length of the horizontal 
error bars is proportional (with a threshold for the 
narrowest ones) to the FWHM of the line, and 
the vertical error bars show the uncertainty in the 
measurement. The two dotted lines are the ratios for 
£(4405 - 5495) and Ay. 


good data. In that sense, the 8621 DIB appears 
to be of cr type i.e. it samples the diffuse, UV- 
exposed ISM but not the dense, UV-shielded 
ISM. 

The two sightlines have different val¬ 
ues of IT(5797)/IT(5780) (0.34 and 0.28 for 
LS III +46 11 and LS III +46 12, respectively) 
but similar values of 7 ^ 5495 , which is in ap¬ 
parent contradiction with the correlation be¬ 
tween 1T(5797)/1T(5780) and 775495 previously 
discussed. However, that correlation is de¬ 
termined mostly by the extremes (sightlines 
with large values of 7^5495 and sightlines with 
large values of VT(5797)/VT(5780)), which is 
the reason why the Pearson coeffcient is not 
very large. For intermediate values Fig. |5] 


shows that it is possible to have variations 
in 1T(5797)/1T(5780) without changing 7 ^ 5495 . 
In other words, the DIB-carrier population 
seems to change continuously with UV-field 
strength/density while for 7^5495 the variation 
seems to happen mostly at the extremes of the 
distribution (in Hit regions and in very dense 
clouds), with an intermediate “typical” range 
where 7^5495 is nearly constant. 


5. Conclusions 

- The M14 family of extinction laws pro¬ 
vides better optical+NIR fits than either 
F99 or CCM for both the Milky Way and 
30 Doradus. 

- 775495 and VT(5797)/VT(5780) are anticorre¬ 
lated, indicating that extreme ^ clouds have 
a larger fraction of small dust grains than 
extreme cr clouds. However, in the inter¬ 
mediate region between the extremes it is 
possible to have clouds with the same grain 
distribution but different proportions of ^ 
and cr DIBs. 

- The 8621 DIB is strongly correlated with 
77(4405 - 5495) up to Ay ~ 6 and the pre¬ 
liminary evidence for some sightlines point 
towards a cr character. 
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Fig. 9. Model for the ISM clouds present in the LS III +46 11 and LS III +46 12 sightlines. Cloud ci is of 
low density, not associated with the cluster, and affects both sightlines similarly. Cloud <T 2 is also of low 
density, is the skin of the cloud associated with the cluster, and is longer along the LS III +46 12 sightline. 
Cloud f is of high density, is the core of the cloud associated with the cluster, and affects LS III +46 11 
exclusively (or, at least, to a much larger degree than LS III +46 12). 
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